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A dissipative quantum system is treated here by coupling 
it with a heat bath of harmonic oscillators. Through quan- 
tum Langevin equations and Ehrenfest's theorem, we estab- 
lish explicitly the quantum Duffing equations with a double- 
well potential chosen. A quantum noise term appears the only 
driving force in dynamics. Numerical studies show that the 
chaotic attractor exists in this system while chaos is certainly 
forbidden in the classical counterpart. 

PACS numbers: 05.45.+b,03.65.Sq,05.40.+j,42.50.Lc 



Quantum chaos of Hamiltonian systems has been stud- 
ied extensively [jO-H. By contrast, very little work has 
been done in looking at quantum chaos of dissipative sys- 
tems. Many quantum mechanical systems (e.g., SQUID 
with Josephson junction, an atom in a cavity of electro- 
magnetic fields, NMR quantum measurements), however, 
are neither isolated nor Hamiltonian. These interact 
with their environment and thus are open in general and 
noisy and dissipative. Dissipation is relatively difficult 
to treat in a quantum system since it seems inevitable to 
deal with stochastic processes via, for instance, the com- 
monly used master equation that is extremely difficult to 
solve numerically. Recently, Spiller and Ralph j|] stud- 
ied a damped and driven non-linear oscillator by using 
the quantum state diffusion model based on the assumed 
master equation of Lindblad form ||. They simulated 
the behavior of one member of the ensemble with a sin- 
gle environment operator and found that the quantum 
noise "kicks" the motion between the chaotic and the 
periodic behavior and so smears out any fractal struc- 
ture. In general, the quantum state diffusion method 
often demands an approximation of the effects of the 
complicated environment by simple operators and always 
has a problem regarding the dimension of the environ- 
ment Brun Q has tried to derive the quantum ver- 
sion of the forced and damped Duffing oscillator through 
the decoherence approach which involves the path inte- 
gral. As he noticed, performing calculations in the low- 
temperature limit is extremely difficult. For the high 
temperature limit, Brun also suggested that the quantum 
noise "smears out" the quantum maps based on Wigner 
distributions and that numerical computation cannot be 
efficiently done because of the necessity of enumerating 
all the possible histories and elements of the decoherence 
functional. 

The dissipative quantum system has been studied by 
means of several different theories, for instance, the in- 



fluence functional approach of Feynman and Vernon ||] 
and its application to Brownian motion , the quantum 
Langevin equations pifl , the master equations |^| f^lTTf . 
etc. |T^ |. Among these, the quantum Langevin formal- 
ism is more interesting for us since the dynamics of quan- 
tum operators for the system can be explicitly given in a 
fashion quite similar to the classical one. Most of above 
theories considered a generic model of the system-bath 
interaction. Recently, Pattanayak and Schieve || pro- 
posed a derivation of semiclassical dynamics directly from 
Heisenberg equations of motion via Ehrenfest's theorem 
where ostensible (apparent) quantum chaos was found 
for the conservative (Hamiltonian) system with a double- 
well potential. In this paper, we will study a dissipative 
quantum model similar to that of Refs. (^[fl[] and write 
explicitly the quantum Duffing equations for a system 
of double-well potential by means of Langevin formal- 
ism and Ehrenfest's theorem. Numerical results are now 
possible to obtain and show that chaotic attractors ro- 
bustly exist in this model whereas chaos is impossible 
for the corresponding classical Duffing system due to the 
absence of an external driving force. 

To do this, let us consider a particle of unit mass mov- 
ing in a one-dimensional time-independent bounded po- 
tential V(Q) with its Hamiltonian, H sys = ip 2 + V(Q). 
To introduce dissipation, we may take the system H sys 
linearly interacting with an external "heat bath" of many 
degrees of freedom, which here is assumed an assembly of 
harmonic oscillators [PHTll . The familiar example of this 
model is a system of an atom interacting with a bath of 
equilibrium photons. Then the complete Hamiltonian is 

to, 
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in which we have assumed that all harmonic oscillators 
possess the same mass m but may have different fre- 
quency u n and that C n is the coupling constant be- 
tween the system and the nth oscillator. The equal 
times commutation relations implicit in (|l|) are [Q, P] = 

ih,[Q,X n ] = [Q,p n ] = [P,%n] = [P,Pn] = 0,[xi,X n ] = 

IPhPn] — and = ihSi n . Now it is straightforward 

to write down the Heisenberg equations of motion for 
the complete system. Then we may apply the standard 
procedure of Ref. [ p"l"| to derive the quantum Langevin 
equations. Assuming a continuous frequency distribution 
g{oj) of harmonic oscillators and using the first Markov 
approximation 
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where 7 is assumed constant, we then have 

> = -HQ)-7^ + e(i), 
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(0) = (P), 

(P) = -(V'(Q)) 
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where £(f) is the quantum noise operator due to the 
heat bath and 7 defined by ^ clearly represents the 
constant damping coefficient. Suppose that the system 
and the bath are initially independent (at f — > —00) so 
that the complete density operator can be written into 
P = Psys <8> Pbi and that the bath is initially thermal, 
Pb ~ exp(— Hb/ksT). Then, £(t) has the properties (see, 
for example, pj) 



We now expand the equations around the centroid by 
using the identity F(Q) = Y, n F {n) Q n /n\, where fW = 
d n F(Q)/dQ n \Q = (Q}, and for a double-well potential of 

V{Q) — —\aQ 2 + jbQ 4 obtain the closed system of the 
stochastic differential equations 

Q = P, 

P = aQ-bQ 3 - SbQhp - 1 P + 77(f), 
fi = a, 

l + a 2 
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(4a^lere we have written Q, P for (Q), (P). We remind the 
fiLO 1 1 \ 1 reader that all physical observables must be obtained by 

UJ0JCO ^2ksT 0SUJ averaging over stochastic noise 77(f). When damping con- 

stant 7 is set to zero, the system is Hamiltonian with 



in which [•••]+ denotes an anticommutator. The average 
(• • •) is over all bath variables. The operator nature of 
£ (f) can be reduced by using the strategy of the adjoint 
commutative representation [13]. We can define a new 
operator 77(f) by ij(t)p(t r ) = ^[£(f), p(t')]+ for all f and 
f', which yields [77(f), 77(f)] = 0. This means that r){t) is 
a c-number function of time. Let us replace the opera- 
tor £(f) by the c-number 77(f) in (||). Its 1- and 2-point 
correlation functions are given by (Q). In the following 
we will restrict ourselves in low temperature limit (not 
discussed in M), T — s- 0, in which quantum effects are 
most important. In this limit the noise 77(f) is given by 
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6hQ 2 p 2 ), for which p — > p 2 and a — » 2p-K in (Q). The 
above equations of motion can be reduced to those of 
Ref. H , in which Hamiltonian semiquantal chaos was re- 
ported. 

Also, in the classical limit Ti — > the first two equations 
of (0) decouple from the fluctuation variables (p, a) and 
the noise term vanishes as T — > [^9|. The well-known 
classical Duffing equations without external driving are 
recovered (for T — > 0) 



(v(t)) = 0, 

{T}(t)rj(t')) = — / dujujcosu{t' - t) 



1 

'T(f-f') 2 ' 
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(5b) 



Q = p, 

P = aQ-bQ 3 - jP. 



(8a) 
(8b) 



To obtain this, an exponential cutoff (e~ euJ ) was used, 
letting e — > after integration E3 ]. 

While the dependence of bath operators has been elim- 
inated in it, the quantum Langevin equations (^) are 
still operator equations and are as difficult as that in 
the deterministic quantum mechanics. To make fur- 
ther progress, we make an additional assumption that 
the wave packet of the system can be described by 
the squeezed coherent state [H. Then we have the 
relations || (Q 2rn ) = {2m)\{hp) m /m\2 rn , (Q 2m+1 ) = 

°j ^ = ri (! + ° 2 )/ 4 ^' (QP + PQ) = ha > where 

0=0 — (O) henceforth and (• • ■) denotes the expectation 
value. It has been shown that these are exactly equiva- 
lent to those derived from the generalized Gaussian wave 
functions ]T^-jl8[|. The equations of motion for the cen- 
troid of a wave packet representing the particle are given 
from (|^) by 



This shows that (0) are indeed the quantum analog of 
Duffing equations with the presence of a quantum noise 
term, which serves as a Langevin driving force. Now a 
question perhaps rises upon how effective our equations 
(0) could be to give the description of exact quantum dy- 
namics. The procedure of using Gaussian wavepacket to 
describe the motion of a particle may be subject to large 
errors and will even break down after a certain time scale 
pOf . This has yet to be decided. It has been argued, 
for example, by Heller |T^] who has extensively studied 
quantum chaos without noise that the generalized Gaus- 
sian approximation works extremely well. Also, recently 
Ashkenazy et al Q have computed the time development 
of the wave function in the presence of a potential bar- 
rier in a bounded well for a long time. They numerically 
confirmed the appearance of Hamiltonian chaos due to 
tunneling first suggested in Ref. § . 

In order to proceed with simulations of the stochas- 
tic differential equations ((?]), we write 77(f) = 

dto v{to)\J^^-e~ %UJt , where v(u>) is a random func- 
tion with the properties: (v(u))) = 0, {v((jS)v{uj')) = 
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5(u> — u>'), v{lo) — v(—ui). One can easily check that 
rj(t) generated in this way satisfies (||). It follows that 
(0) can be numerically solved for each realization of the 
random process ry(t) by the Runge-Kutta method for the 
ordinary differential equations once the random sequence 
of v (u>) has been generated. Quantum chaotic attractors 
are then found for weak damping after some evolution 
time (typically around 1000 in our simulations). For one 
realization of stochastic process r](t), Fig. |l](a) shows the 
structure of the quantum chaotic attractor in the phase 
space. The fact that the attractor diffuses out due to 
the quantum noise agrees with the previous results pH . 
In particular, one may compare this attractor with that 
of Ref. § for the case of strong noise. Lyapunov expo- 
nents As and the fractal dimension D f are calculated for 
each realization by using the standard method as given 
by Ref. [^TJ and, as seen in Fig. 1(b), are saturated af- 
ter time of the order t — 10000 with a small residual 
oscillation less than 1% (Xiargest — 0.124500). For this 
case, we furthur checked the largest Lyapunov exponent 
by using an alternative computation algorithm [p2[ and 
found that the Lyapunov exponents Xiargest — 0.124470 
agreeing with the above number up to an accuracy of 1%. 

We have computed 1000 samples (realizations) syn- 
chronously and then have obtained the distribution of 
Lyapunov exponents A and a probability map (see Fig. ||). 
As seen in Fig. |](a), the largest Lyapunov exponents for 
all realizations are conclusively positive. Note that the 
fractal dimension Df is very close to four due to the 
weak damping (7 = 0.002) since £>/ 4 as 7 — > where 
the system becomes Hamiltonian. Now let V(Q,t) de- 
note the probability of the system at position Q at time 
t. We make use of a probability map j2^| defined by 
V(Q,t) — > V(Q,t + t) for constant Q and r to further 
characterize the behavior of the noisy quantum system. 
It has been shown [^| that for regular behavior the prob- 
ability function V(Q,t) does not depend on time once 
the motion of the system is stable and therefore the map 
V(Q,t) — » V(Q,t + t) should only consist of a single 
point. Hence, Fig. ||(b) justifies the chaotic behavior of 
this noisy quantum system as well p4j . By contrast, it is 
well known that for the corresponding classical Duffing 
equations flq) only point attractors, describing regular 
motions, can be allowed due to the absence of the ex- 
ternal driving force. While it has been suggested ||^|] 
that the tunneling effect induces chaotic behavior in a 
Hamiltonian quantum system, the dynamical effect of the 
quantum noise in (0) should have played a crucial role to 
form the stable chaotic attractor in our dissipative quan- 
tum system. Otherwise, the system has to be attracted to 
the bottom of either well with zero momentum. In other 
words, we report here that it is the quantum noise that 
leads to the chaotic attractor. Indeed, this is somewhat 
reminiscent of the classical fact that either multiplicative 
or additive noise may induce homoclinic crossing and so 
chaos, as suggested by Schieve, Bulsara and Jacobs pj| 



in their studies for classical stochastic chaos. 

To conclude we note that the quantum Duffing equa- 
tions in low temperature limit have been explicitly es- 
tablished from the quantum Langevin equations. These 
equations manifestly display great advantage of numer- 
ical computation in comparison with others Nu- 
merical results show that these equations exhibit the sta- 
ble chaotic attractor for weak damping while as already 
known the classical counterparts certainly forbid chaos 
due to the absence of an external driving force. To our 
knowledge, this is the first study of the dynamical behav- 
ior of the dissipative quantum system without external 
driving showing the quantum chaotic attractor for such 
a system. More detailed studies shall be presented else- 
where [^9). 

We should bear in mind that some assumptions have 
been made. First, all results are based on a simplified 
theoretical system-plus-environment model for which we 
assumed that: i) the heat bath consists of an assembly 
of harmonic oscillators; ii) there is a continuous distri- 
bution of oscillator frequencies; Hi) the coupling of the 
system to the bath operators is linear and the coupling 
constant is a smooth function of oscillator frequency; 
and iv) the stochastic process is Markovian. All these 
are quite well-known and usual in the study of a dissi- 
pative quantum system (see, for example, Refs. [p|-[il"|| ) . 
Second, the squeezed coherent state (or equivalently the 
generalized Gaussian wave-packet) has been used to ap- 
proximate the true wave function of the system. The full 
quantum phase space is thus restricted into a truncated 
"semiquantal" phase space p|,p0| . It has been shown by 
Ashkenazy et al. [0 via computer simulation that this ap- 
proximation does not break down for a long time. Fully 
understanding its validity is still an open task. 

The authors would like to thank Arjendu K. Pat- 
tanayak for many stimulating discussion during his re- 
cent visit. 
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FIG. 1. Quantum chaotic attractor for one realization with 
parameters: h = 1, a = 10, 6 = 4, 7 = 0.002. (a) The 
Poincare section Q-P taken by a — ± 0.0015 and a > 0; 
(b) The time evolution of the largest Lyapunov exponent A. 
For this realization, the sum of all Lyapunov exponents in 
phase space ^\ A,: ~ —4.011 x 10 -3 and the fractal dimension 
D f ~ 3.968. 
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FIG. 2. The simulation for 1000 realizations with the same 
parameters as in Fig. |l|. (a) The distribution D(X) of the 
largest Lyapunov exponents with J d\D(\) — 1. The mean 
Aiargost ^ 0.127, D f ~ 3.969 and J2i X * - -3.957xl0" 3 . Lya- 
punov exponents are calculated by means of Wolf et al pl| . 
(b) The probability map, taking constant Q = 1.581 ± 0.05 
(near the bottom of one potential well) and r = 5. This map 
starts after a relaxation time of 10000. 
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